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SUMMARY 


A linear, inviscid, subsonic compressible flow theory is formulated for 
predicting the aerodynamic characteristics of upper-surface-blowing 
configurations. The effect of the thick jet is represented by a two-vortex- 
sheet model in order to account for the Mach number nonuniformity. The wing 
loading with the jet interaction effects is computed by satisfying boundary 
conditions on the wing and the jet surfaces. The vortex model is discussed 
in detail. 


INTRODUCTION 


In upper-surface-blowing (udB) configurations, the low-pressure-ratio 
jet from high by-pass ratio turbofan engines blows directly on the wing 
upper surface. As the jet is relatively thick, being of the order of 10% 
of the local chord, the conventional thin jet flap theory has been found to 
be inadequate to predict the high lift (ref. 1). This means that additional 
lift may come from the interaction between the wing flow and the thick jet 
which has higher dynamic pressure than the freestream. Of course, the jet 
entrainment will increase the lift also, mainly through producing the Coanda 
effect. t 

In this paper, a theoretical method will be described for predicting 
the interaction effects due to nonuniformity in Mach numbers and dynamic 
pressures in the flow field. The inviscid, linear, subsonic compressible 
flow theory is assumed. 


SYMBOLS 


total lift coefficient (circulation lift plus jet reaction lift) 


* This work was supported by NASA Langley Grant NSG 1139. 
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ACl difference in lift coefficients with jet on and off 
induced circulation lift coefficient 
Cjj jet momentum coefficient 

c local chord length, m(ft) 

Cj length of trailing jet included in the analysis, m(ft) 

M number of integration points or Mach number 

N number of chrodwlse integration points 

Nj number of streamwise vortices on the trailing jet 
n unit vector normal to jet surface 

R radius of curvature, m(ft) 

s coordinate tangential to the jet surface 

T = Po/ Pj 

tj jet thickness, m(ft) 

u nondlmenslonal backwash 

V velocity, m/sec (ft/sec) 

x,y,z rectangular coordinates, with x positive downstream, y positive spanwise- 
to the right and z positive upward 

a angle of attack, deg. 

V nondimensional vortex density 

6f flap angle, deg. 

6j jet-deflection angle, deg. 

V - V„/Vj 

p' * y cos a 

p density, kg/m^ (slugs/ft3) 

(fi velocity potential, m^/8ec(ft^/sec) 

ij» nondimensional additional perturbation velocity potential 
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Subscripts: 


c chordwise 

j jet flow 

o outer flow 

s spanwise 

w wing 

°° freestreain 

METHOD OF ANALYSIS 
Formulation of the Problem 

The perturbation flow fields Inside and outside the jet are assumed to 
be governed by the Prandti-Glauert equations with Mj and Mq, respectively. 
The solutions of these equations must satisfy the boundary conditions that 
the jet surface is a stream surface and the static pressure must be con- 
tinuous across it, in addition to the usual wing tangency condition. In the 
linear theory, these conditions can be written as (ref.l) 


34>o 

3ij>j 

-V„.n(l-p') 

(jet stream surface condition) 

(1) 

8n 

3n 

io.t 



aij 

ds 

- - T(p’ 

3(i) 

9 “O 

■ 0 

(jet static pressure continuity) 

(2) 

9^0 

dz 

3^c 

■ “ 3x 

- tan OL 

(wing tangency) 

(3) 

where 





T « 

Po/Pj 



(4) 


and (pQ and (j)j are related to the dimensional velocity potentials (p^ avj (pj as 

Po * 4>0 Voo cos a, pj = pj Vj C5) 

Since the above problem is linear, it can be decomposed into a wing-alone 
case with potential and the interaction case with additional potential i|^. 
Let _ _ 

<l>o = Ko (Mo) + (Mo) (6) 

= iwj (Mj) + 'Pj (Mj) (7) 

where 4>wo *^wj satisfy the Prandti-Glauert equations with respective 
Mach numbers and the following boundary conditions 
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tan a 


( 8 ) 


L 


3't'wo 

d Zp 

(Mo) ^ 

3z 

3x 


3zc 

Bz 

(Mj) . ^ 


tan a 


(9) 


Substitution of equations (6)-(7) into equations (l)-(3) gives 
Bipo -nd-y') 


9'J^i 

^ -T(y-) 


2 


S'I'o 

Bs 


3<f>wj 

Bs 


+ T(y')2 


3<t>wo 

3s 


3<l^o 


3z 


0 


( 10 ) 

( 11 ) 

( 12 ) 


The above equations indicate that there are jumps 
tangential velocities across the jet surface. If 
Mq = Mj and T * 1) , these equations show that the 
potentials will be automatically zero. 


in normal velocities and 
there is no jet (i.e., p*“l, 
additional perturbation 


Vortex Model 


In order to satisfy equations (10)“(11) simultaneously, two vortex 
sheets are introduced on the jet surface. These vortex distributions will in-- 
duce normal and tangential velocities in the flew field. To evaluate these 
induced velocities on the boundaries so that boundary conditions, equations 
(10)-(12) , can be satisfied, the ‘‘nduced v( 'ocity integrals are reduced to 
finite sums through the Quasi Vortex-Lattice Method (Quasi VLM) (ref. 2). 

The resulting discretized vortex arrangement is as shown in fig. 1. Note that 
the wing vortices directly below the jet surface are arranged so that they 
coincide with the jet vortices in location. Furthermore, the shaded region 
represents part of the nacelle. Since the present computer program does not 
include the nacelle model, special care must be exercised in this region. If 
vortex representation of the jet starts from the exit, instead of the leading 
edge, it has been found that the jet would Induce a leading-edge type loading 
in the middle of the wing. Therefore, the vortex distribution must be extend- 
ed to the wing leading edge, but with freestream conditions assumed in the 
nacelle region in the computation. 

With the vortex arrangement made, the required induced velocities can be 
computed and substituted into equations (10)-(12) for the solution of unknown 
vortex strengths. These equations are satisfied only at a discrete number of 
control points which are chosen according to the *’seml-circle method** as de- 
scribed in the Quasi VLM (ref. 2). These calculations are mostly straightfor- 
ward, except in computing the induced tangential velocities (i.e., backwash) 
and representing the jet flap effect. These will be discussed below. 
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Computation of Induced Tangential Velocities 


The induced tangential velocities on the jet surface are needed to sat- 
isfy equation (11). At any control point, the induced tangential velocity di 
to the jet vortex distribution in its neighborhood is simply equal to the 
vortex density at that control point. Since the vortex density at the control 
point does not appear in the formulation, it is necessary to express it in 
terms of those at the vortex points through interpolation, such as Lagrangian 
interpolation. However, the contribution from those vortices not on the same 
plane with the control point can not be accurately computed in the usual man- 
ner as in computing the induced normal velocities. To illustrate, consider the 
backwash expression in the two-dimensional flow: 


u(x,z) 


_Z ri Y(x_*_)dx' 

2ir J Q (x-x')^ + z 2 


(13) 


As z 0, the integrand of equation (13) will have a second-order singularity. 
The usual method is not accurate in treating such singularity. Therefore, 
equation (13) should be rewritten as 


u(x,z) 


rl Y(x») - Y(x) 

2tt J 0 (x-x')2 + z2 


dx' + 


y(x)z r ^ 
2tt Jo 


dx' 

(x-x')^ + 2 ^ 



£ 

N 



Y(6k)-Y(x) 
(x-xk)2 + z2 


sin 0k + 


2Y(x) £ 
Air M 


M sin 6j 


7 


(lAa) 


where 


Xk = (1 - cos 9k) , 0k 

Xj = (1 - cos 6j), 


2k - 1 
2N 

2j - 1 
2M 





(lAb) 


and 

M = 2^N 


(lAc) 


for interdigitation between the control and integration points in the last 
summation and p is an arbitrary integer greater than 1. 


In the implementation of equation (lAa) in the present computer program, 
it is assumed that M*8N if N>6 and M“16N if 6. To compute the tangential 
velocity at a jCt control point due to wing vortices, equation (lAa) is ap- 
plied to three wing vortex strips in the immediate neighborhood of the control 
point. This is indicated in fig. 2. On the other hand, if the control point 
is on the jet upper surface, the effect due to all vortices on the jet lower 
surface in the same streamwlse section is computed with equation (lAa). 

Similar principle is applicable if the control point is on the lower surface. 
This is Illustrated in fig. 3. The accuracy of equation (lAa) has been 
illustrated in ref. 3. 
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Incorporation of Jet Flap Effect in Thick Jet Tneory 


Due to Coanda effect, the USB jet will follow the wing surface and leave 
the wing trailing edge at an angle relative to the ch^rd t^ produce the jet 
flap effect. This is to produce varying but unknown •• n in equation (10) 
along the trailing jet. ^ 

Since 



dz 

dx 


+ tan a 


(15) 


it is necessary to relate to the unknown jet vortex density. As shown in 
fig* irrotationallty of the jet flow implies that 




B(pj 


or 


R 


T(p*)2[(— 


3«l»o 3'^wo 




r 


2 . - 

R dx^ (17) 

where the subscripts 1 and 2 denote the upper and the lower jet surfaces « 
respectively, at the section under con&iderat:Ion. Equations (16)-(17) give 
the following initial value problem for determining ^ : 

dx 




z 

z(t.e.) 


dx^*^*®'^ 


f(x,y) 

0 



(18) 


To integrate the above problem for , a finite jet length Cj is first chosen 
(It is not necessary tc include infinite length of jet in the numerical 
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method as far as the wing loading is concerned) . Then let 

X ■ JCt.e. + -^ (1 “ cos e) 

It follows that 


and after integration once. 



sin 0f(b,y)d0 


(19) 


Ci p i“l . 

« -tj 6j + (A0) ^ I sin 0kf(0k.y) + ^in 6if(0i,y) (20) 

dz k X 

by trapezoidal rule. With determined by equation (20) in terms of 
f(B,y) and hence, — 2. through equation (16), it is possible to incorporate 

dS 

these terms with unknown vortex densities to the left hand side of 
equation (10) before equations (10)-(12) are solved. 


Note that equations (10)-(12) are solved by the vector method of Purcell 
(ref. 4) which processes row by row of the augmented matrix in solving the 
equations. Since the tangential velocities are needed in equation (10) with 
the jet flap effect (see equation 16), equation (11) must be processed first 
with the tangential velocities computed there stored on file before 
equation (10) can be solved. 


SOME NUMERICAL RESULTS 


Before the method can be applied to any configurations, it is important 
to know how the discretized vortices should be arranged to produce reliable 
results. In the following, some convergence study with respect to the vortex 
arrangement for the configuration used by Phelps, et al. (ref. 5) will be pre- 
seated. Fig. 5 shows the effects of number cf vortices on the trailing jet 
with Cj ■ Ic and 2c. It is seen that 6 vortices (Nj * 6) with Cj = Ic appear 
to be sufficient in appll'^ations. If Cj is increased to 2c, more vortices 
would be needed to provide convergence. Fig. 6 indicates the rapid convergence 
of ACl^ with respect to the number of wing spanwise vortex strips. 

The above method is now applied to the configuration of ref. 5 with 30® 
full span flap. The jet-deflection angle is taken to be the angle of flap 
extension relative to the chord line. It is found to be 4 6.7®. The Cj^ values 
were computed by adding predicted AC^ to the experimef»tal jet-off C|^. The re- 
sults are «?hovm in fig. 7. The predicted values are slightly higher possibly 
because the estimated jet angle was too high. On the other hand, the thin jet 
flap theory by the present method (ref. 1) v;ould underestimate the lift. 
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To illustrate the importance of improved tangential velocity evaluation 
(see equation 14a), the configuration of ref. 6 with 6f = 20^ and 40^ is 
used in fig. 8. It is seen that without the improved integration technique, 
the predicted lift would be too high. 


CONCLUDING REMARKS 


A vortex model for the USB jet-wing interaction has been described. The 
model consists of using two vortex sheets on the jet surface to account for 
Mach number nonunifonnity and differences in jet and freestream dynamic 
pressures. The rate of numerical convergence with respect to the number of 
vortices used appeared to be reasonably rapid. Comparison of the predicted 
results with some available data showed much better agreement than the thin 
jet flap theory. 
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Figure 1.- Vortex model for computing jet-interaction effects. 



Figure 2.- Region (shaded area) of wing vortex strips subject to improved 
integration procedure of equation (14a) for computing u-velocity at a 
jet control point. 




Figure 3.- Region (shaded area) of jet vortex strips on the lower surface 
subject to improved integration procedure of equation (14a) for com- 
puting u-velocity at a jet control point. 



Figure 4.- Geometry of a jet cross section 
for formulating jet flow Irrotationality 
condition. 
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Figure 7 


O EXPERIMENT (REF. 5) 
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- Comparison of predicted lift curves with experimental data of 

ref. 4. C = 0.8 X 2.095. 

\i 



Figure 8.- Effect of backwash evaluation on predicted lift. C,. = 2. 



